Hardness of T-carbon: Density functional theory calculations 
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We revisit and interpret the mechanical properties of the recently proposed allotrope of carbon, 
T-carbon [Sheng et at, Phys. Rev. Lett., 106, 155703 (2011)], using density functional theory in 
combination with different empirical hardness models. In contrast with the early estimation based 
on the Gao's model, which attributes to T-carbon an high Vickers hardness of 61 GPa comparable 
to that of superhard cubic boron nitride (c-BN), we find that T-carbon is not a superhard mate- 
rial, since its Vickers hardenss does not exceed 10 GPa. Besides providing clear evidence for the 
absence of superhardenss in T-carbon, we discuss the physical reasons behind the failure of Gao's 
and Simiinek and VackM's (SV) models in predicting the hardness of T-carbon, residing on their 
improper treatment of the highly anisotropic distribution of quasi- sp^-like C-C hybrids. A possible 
remedy to the Gao and SV models based on the concept of superatom is suggest, which indeed 
yields a Vickers hardness of about 8 GPa. 

PACS numbers; 64.60.My, 64.70.K-, 62.25.-g, 62.20.Qp 
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Recently, on the basis of first-principles calculations 
Sheng et al. proposed a carbon allotrope which they 
named T-carbon^. Strictly speaking, its actual stability 
needs a highly large negative pressure which is far be- 
yond currently available technologies. Structurally, this 
phase can be obtained by substituting each carbon atom 
in diamond with a carbon tetrahedron (Fig. [T]), and thus 
crystallizes in the same cubic structure of diamond (space 
group Fd3m)) with the carbon atoms at the WyckofF site 
32e (0.0706, 0.0706, 0.0706). It has been noted that T- 
carbon has a large lattice constant of 7.52 A and a low 
bulk modulus of B =169 GPa, only 36.4% of the bulk 
modulus of diamond^. In particular, its equilibrium den- 
sity, 1.50 g/cm'^, is the smallest among diamond (cubic 
and hexagonal diamond}^, graphite^, M-carbon^, bct- 
C^^, H^-carbon^, chiral-carbon^ as well as the newly pro- 
posed dense hpi-, tI12- and iP12-carbon^ phases. This 
results in an highly porous structural pattern, which can 
be viewed as a diamond-like array of superatoms (tetra- 
hedral C4 clusters), as depicted in Fig. [TJ Given this pe- 
culiar clusterized arrangement of atoms exhibiting a quite 
low shear modulus of G = 70 GPa^, it is very surprising 
that T-carbon was predicted to be superhard, with an 
exceptionally high Vickers hardness (H^) of 61.1 GPa^, 
comparable to that of superhard cubic boron nitride (c- 
BN). 

The aim of our present study is to elucidate the ori- 
gin of this anomalous hardness. We do this by exploring 
in details the mechanical properties of T-carbon through 
the application of several different empirical approaches: 
the Gao's formula^, the SV model® and our recently pro- 
posed empirical treatment based on the Pugh's modu- 
lus ratioS.. Our systematic analysis provides an unam- 
biguous and physically sound results: T-carbon is not 
hard. We will show that the conventional application of 
Gao and SV models leads to a much too high Vickers 



hardness, H,^^°=61.1 GPa and GPa, substan- 

tially overestimated with respect to the value obtained 
usmg our formahsm, H^h™=5.6 GPa. The prediction 
of a low Vickers hardness in T-carbon is consistent with 
the estimation of a low shear strength (7.3 GPa along 
the (100)<001> slip system), which represents the up- 
per bound of the mechanical strength. 




FIG. 1: Lattice structure of T-carbon (space group Fd3m). 
By considering each carbon tetrahedron (C4 unit) as an ar- 
tificial superatom, the corresponding structure is isotypic to 
that of diamond. The local environment of each superatom is 
illustrated in the right panel. 

The improper assignments derived by a conventional 
application of the Gao and SV models can be attributed 
to the fact that these two models assume that the chem- 
ical bonds, which are significant for hardness, are dis- 
tributed uniformly in the lattice. But in T-carbon, as 
already pointed out be Sheng et al} , though the car- 
bons atoms are tetrahedrally coordinated and apparently 
resembling a three-dimensional quasi- sp^-like hybrid^, 
their bonds are ordered in an extremely anisotropic and 
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porous framework, highly different from the bonding dis- 
tribution in ideal sp^-hybrid. We propose a remedy to 
cure the limitations of Gao and SV models in dealing with 
anisotropic and porous systems by assuming each carbon 
tetrahedron cluster as an artificial superatom. Indeed, 
this cluster-like approach leads to low Vickers hardness 
in the range of 7-8 GPa, in agreement with the estimated 
value of 5.6 GPa using our proposed model^. 

All calculations were performed using the Vienna ab 
initio Simulation Package (VASP)^ in the framework 
of density functional theory (DFT), and we adopted the 
Perdew, Burke and Ernzerhof approximatiorti^ to treat 
the exchange-correlation kernel. Well converged results 
were obtained using an energy cut-off of 500 eV and a 
k-point grid 11x11x11—. The DFT results were then 
employed as input for the three different hardness empir- 
ical models, with which we have computed the Vickers 
hardness H.^: 

(a) Gao's model^: 

i/Gao ^ 350[(iV2/3)e-1.191/./(rf2.5)] 

where Ne is the electron density of valence electrons per 
A'^, d is the bond length and fi is the ionicity of the 
chemical bond in a crystal scaled by Phillips. As already 
mentioned, this model gives H^'^°=61.1 GPa (Ref. l). 

(b) SV model«: 

(J 

Hy"^ = ■^V^^/{d^^n^^) (2) 

where C is the constant of 1550 and ft is the equilib- 
rium volume of T-carbon. — Zi / Ri represents the ref- 
erence energy, with Zi indicating the valence number of 
element i. For carbon — 4.121 (taken from Ref. ^. nu 
and da are the number of bonds and bonding lengths be- 
tween atom i. In T-carbon, each carbon has four nearest- 
neighbors with two different bonding lengths: three in- 
tratetrahedron carbon-carbon bonding length of 1.502 A 
and one intertetrahedron bonding length of 1.417 A^. 
By using the average bonding length of 1.48075 A we ob- 
tained Hy = (1550/26. 5785)x4.121/(1.48075x4) = 40.5 
GPa, which is 33.5% smaller than the corresponding 
Gao's value. ^ 

(c) Chen's model^. This is the empirical formula which 
we have recently proposed, based on the Pugh's modulus 
ratio k=G/B^: 

This model not only reproduced well the experimental 
values of Vickers hardness of a series of hard materials 
including all experimentally verified superhard materials 
(see Fig. [5] and Table |l|), but also provides a theoret- 
ical foundation of Teter's empirical correlatioiJ^ in its 
simplified form^. 

Before discussing the results for T-carbon we start by 
presenting some general considerations regarding the cal- 
culation of the Vickers hardness and the trustability of 



our proposed model^ . Hardness is a highly complex prop- 
erty, which depends on the loading force and on the qual- 
ity of samples (i.e., presence of defects such as vacancies 
and dislocations). Because Vickers hardness is exper- 
imentally measured as a function of the applied load- 
ing forces, the saturated hardness value (or experimen- 
tal load-invariant indentation hardness) is usually consid- 
ered to be the hardness value of a given material. There- 
fore, the theoretically estimated Vickers hardness within 
Gao's, SV's and Chen's models should be directly com- 
pared to the experimentally saturated hardness value of 
polycrystalline materials. An overview on the experimen- 
tal and theoretical values of H„ for the experimentally 
verified superhard materials (diamond, C-BC2N, c-BN, 
C-BC5, 7-B28) is summarized in Fig. [5] and Tab. HI The 
experimental results are highly scattered, reflecting the 
inherent difficulties in achieving a trustable and precise 
estimation of hardness. For instance, the reported values 
for the hardness of diamond, the archetype superhard 
materials, range from 60 GPa to 120 GPaii^"— . Similar 
trends have been observed for the other two well-known 
superhard materials C-BC2N and c-BN. The most typ- 
ical case is probably ReB2, whose actual hardness has 
been extensively debated32r^ after the first value of its 
Vickers hardness (48 ± 5 GPa at the loading force of 0.49 
N) was reported^i. Depending on different samples, syn- 
thetic methods and measurement technique, the obtained 
values range from 18 to 48 GPa (Table IJ). In contrast to 
experiment, theoretical estimations of the Vickers hard- 
ness given by different models^^^ agree within few GPa, 
including the data obtained by our proposed model (Eq. 
Is}. Overall, the comparative trend displayed in Fig. [2] 
provides robust evidence for the reliability of our pro- 
posed formalisrr^. 

Now, let's turn the attention to T-carbon. By using 
the values of the shear and bulk moduli from Ref 1 as in- 
put {B = 169 GPa and G = 70 GPa) for Eq. |3]we obtain 
a Vickers hardness of 5.6 GPa, dramatically smaller than 
the corresponding Gao (61.1 GPa) and SV (40.5 GPa) es- 
timations. Furthermore, we noted that Sneddon defined 
the concept of ideal elastic hardness by Hid — 2^1™* t) 
where E is Young's modulus, v is Poisson's ratio and 
coi(j) w 0.5 for the standard pyramid indentation"*"^ and 
suggested that the real hardness would be (0.01^0.2)H,irf 
at high loads^. Utilizing this definition and the derived 
E = 185 GPa, the real hardness for T-carbon should be 
in the range from 0.5 GPa to 10 GPa, in agreement with 
our obtained value. In particular, it still needs to note 
that the occurrence of this serious discrepancy among the 
three different methods (Gao's, SV's and Chen's models), 
which is not observed for the other test cases of Fig. [2] 
and Table HI urges for a clarification aiming to discern 
which method provides the more reliable description of 
the hardness of T-carbon and, consequently, to help us 
to answer a naturally arising question: is T-carbon a real 
superhard material? 

An useful concept for understanding strong mechani- 
cal strength - but still relying on elastic properties - is 
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FIG. 2: Vickers hardness as a function of a product (k^G) 
of the squared Pugh's modulus ratio {k = G/B) and shear 
modulus (G). The curve corresponds to the empirical relation 
of Eq. [3] (For other data and more details, see Ref. [^. Elas- 
tic moduli and experimental Vickers hardness are collected in 
Table |T] Note the huge discrepancies among the three theo- 
retical estimations for T-carbon. 



based on ideal shear and tensile strengths^, at which 
a material is getting unstable under direction-dependent 
deformation strains^. To shed some light on the na- 
ture of T-carbon we have thus investigated ideal tensile 
strength along the <001> direction and shear strength 
along the (100)<001> slip system. We found that a ten- 
sile strength of 40.1 GPa along the <001> direction and a 
shear strength of 7.3 GPa in the (100)<001> slip system 
(see Fig. [3]). Therefore, we can conclude that the failure 
mode in T-carbon is dominated by the shear deforma- 
tion type in the (100)<001> slip system. The calculated 
shear stress of 7.3 GPa basically sets the upper bound on 
its mechanical strength at zero pressure^^i^, because the 
ideal strength is the stress where a defect-free crystal be- 
comes unstable and undergoes spontaneous plastic defor- 
mation. It is well-known that the measurement of hard- 
ness has to first encounter the elastic deformation and 
then experience permanent plastic deformation. There- 
fore, it can be conjectured that the hardness of T-carbon 
should not exceed 7.3 GPa. These arguments provide a 
strong support for our estimated Vickers hardness of 5.6 
GPa on the basis of Eq. [S] 

In order to gain further insights on this intricate sub- 
ject and to reach a consistent and satisfactory conclusion 
on the hardness of T-carbon we consider now the rela- 
tion between hardness and brittleness on the basis of the 
Pugh's modulus ratioi^. There is no doubt that all exper- 
imentally verified superhard materials, such as diamond, 
c-BN, C-BC2N, 7-B28 and C-BC5 are intrinsically brittle. 
As shown in Table U the Pugh's modulus ratic^'^ of these 
superhard materials [k = 1.211-1.178 (diamond), 0.999- 



TABLE L Comparison between measured (H^^^) and the- 
oretically computed values of the Vickers' values (in GPa), 
along with available bulk modulus (B, GPa), shear modulus 
(G, GPa) and Pugh's modulus ratio fc= G/B. 
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FIG. 3: DFT calculated ideal tensile and shear strengths of 
T-carbon. 



1.053 (c-BN), 1.107-1.091 (c-BCsN), 1.054 (7-623) and 
1.048 (c-BCs)) are larger than 1.0. They clearly obey to 
the empirical relation that considers the Pugh's modulus 
ratio as an indicator of the brittleness or ductility of ma- 
terials. The higher k the more brittle (and less ductile) 
the material is. Pugh still proposed, when k is larger 
than 0.571 the materials are brittle and with k being less 
than 0.571 the materials are ductile ^2.. This relation has 
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been extensively applied not only to metals and alloys 
but also to high-strength materials. In the case of T- 
carbon, the calculated Pugh's modulus ratio k = 0.414 is 
smaller than 0.571, clearly in the range of ductility. The 
ductile behavior of T-carbon is a further indication of its 
non-superhardness. 

On the basis of the above consideration we can now 
understand why T-carbon is not an superhard material. 
One common feature of superhard materials is that they 
not only need a three-dimensional network composed of 
short, strong, and covalent bonds^ but also have a uni- 
form distribution of strong covalent bonds. The proto- 
typical example is diamond, which is characterized by 
an isotropic array of tetrahedrally bonded sp^^ carbon 
atoms. Conversely, in soft graphite the sp^-type covalent 
bonds, though strong, are localized in two-dimensional 
sheets. At first glance, T-carbon seems to be a good 
candidate for superhardness since each carbon atom has 
four nearest-neighboring carbon tetrahedrally bonded by 
short and strong carbon-carbon covalent bonds. How- 
ever, due to the extreme-anisotropic arrangement of these 
carbon-carbon bonds and the associated formation of a 
large proportion of porosity in lattice space as well as the 
low density of bonds, the framework of T-carbon be more 
easily bendable in comparison with that of diamond, as 
manifested by its low shear strength. 

Having this in mind, we can look back Gao's and SV 
models. Although these two models perform very well 
for many hard materials, they deliver questionable num- 
bers for T-carbon in sharp contrast with our findings, 
as we have documented above. The reason for this ap- 
parent failure is that in these two models all bonds are 
treated as uniformally distributed in the lattice space. 
Clearly, this constrain will not affect the predictions for 
isotropic material but it will be inadequate to describe 
the hardness of extremely-anisotropic compounds such 
as T-carbon. However, if we give a closer look to each 
individual C4 tetrahedron unit (see Fig. 1), the distribu- 
tion of six strong carbon-carbon covalent bonds within 
each C4 unit is highly dense. It is therefore trustfully 
expected, that the Vickers hardness of each individual 
C4 unit can be comparable (or even harder) to that of 
diamond because its bonds density and strengths within 
each C4 unit are higher than those of diamond. The 
strength and rigidity of each individual C4 unit appear 
to be such strong that it cannot be broken easily. Based 
on this fact, in order to render Gao's and SV's meth- 



ods applicable to T-carbon, each carbon tetrahedron (C4 
unit) is considered to be an artificial superatom (See right 
panel of Fig. [T]). The cubic unit cell of T-carbon consists 
of eight superatoms and each superatom has four nearest 
neighbors with the bonding length of d — 3.257 A . In 
terms of Gao's and SV's methods, this distance d should 
be the bonding length between exact atomic positions 
with positively charge cores, representing the real force 
center of each atom. Based on our assumption, the d 
distance is defined as the spatial separation between two 
nearest neighbor superatom positions, d— 3.257 A. Al- 
though it remains disputable whether the center of mass 
of the C4 superatom could be assigned its real force cen- 
ter (thus allowing the applicability of Gao's and SV's 
models) the high strength and rigidity of each individ- 
ual C4 unit manifested by the dense and strong carbon- 
carbon bonds seem to validate this assumption of d dis- 
tance. Obviously, each superatom contains 16 valence 
electrons, and Ne = 8/26.61 = 0.3. By inserting these 
values of d and Ne in Gao's formula (Eq. [T]), we derive 
a Vickers hardness of 8.2 GPa, in agreement with our 
value of 5.6 GPa. To apply the same adjustment to the 
SV model one needs to define the crucial parameter Ri for 
the superatom. From our first-principles calculations, it 
can be inferred that Ri = 2.32 A represents the optimum 
radius containing all 16 valence electrons for each super- 
atom. By inserting = 16/2.32 = 6.896 in Eq. Ha Vick- 
ers hardness of 7.7 GPa is obtained, again in agreement 
with our analysis. Within this superatom approach, all 
three methods discussed in the present paper convey the 
same answer: T-carbon is not superhard. The anomalous 
behavior of Gao and SV models observed in Figl2] for T- 
carbon is cured and the general agreement among the 
three Gao, SV and Chen models is re-established. This 
provides clear evidence that the hardness of T-carbon 
should not exceed 10 GPa. 
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